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The technique, referred to as simple enhancement, is a 
special case of the Gandin (1963) optimum interpolation 
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observation. Trial values are obtained from a climatology 
and a synoptic analysis/forecast system provided by the Fleet 
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is calculated by adding an observed anomaly (i.e., observation 
minus trial value) to the trial value at the desired location. 

Calculations of mean and RMS errors indicate that simple 
enhancement can provide a closer estimate to actual conditions 
than unenhanced climatology. The mixed layer depth cannot be 
extrapolated accurately to new locations presumably due to 
mesoscale eddies, fronts, internal waves and small scale fluc- 
tuations at the base of the mixed layer. Also, the choice of 
the trial value used is not critical. Experiments at different 
locations and seasons 'would be required for a complete assess- 
ment of the application to ASW operations. 
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I. INTRODUCTION 



Accurate ocean thermal structure information is an impor- 
tant requirement necessary to support U.S. Navy anti-submarine 
warfare (ASW) operations. Estimates of oceanic conditions 
that are to be used to generate sound velocity profiles and 
acoustic range predictions can be obtained from climatology or 
objective analyses. These sources of information are usually 
represented digitally on coarse grids which limit their 
spatial resolution. Neither climatologies nor objective 
analyses can depict the true oceanic conditions. 

The best possible source of local information for a ship 
at sea is an observation. The oceanic conditions at a remote 
point, however, are often essential for mission planning. 

The question that arises at sea is how can a local observa- 
tion be used to estimate conditions at some distant point. 

In this thesis, a simple technique for accomplishing this 
will be developed and tested. 

A. EXTRAPOLATING A LOCAL OBSERVATION TO A REMOTE POINT 

Some attempts have been made operationally to "enhance" 
a climatological temperature profile at a remote point by 
adding the departure from climatology determined from a 
local observation to the climatology at a remote point. 

Such a technique has been reported to provide useful esti- 
mates for computing acoustic ranges in Navy fleet exercises 
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(LT D. Pedneau, personal communication) . This type of 
extrapolation is also used in the Expanded Ocean Thermal 
Structure analysis (EOTS) (Holl et al., 1979) run daily at 
Fleet Numerical Oceanography Center (FNOC) . A preliminary 
step in EOTS is to estimate values at grid points from near- 
by observations using this method. 

Various methods of extrapolation also constitute the 
fundamental approach used in other objective analyses. 

Simple methods such as used by Druyan (1972) involve assigning 
values to grid points by weighting the observation inversely 
proportional to the square of the distance from the grid 
point. Other more sophisticated methods require the use of 
a "trial value" at any arbitrary location (Kruger, 1969). 

The concept of a "trial value" is used in the same context 
that a climatology is used in this thesis. In these applica- 
tions, the deviation of the trial value from the observation 
(climatological anomaly as used here) is multiplied by a 
weighting factor before adding to the trial value (clima- 
tology) at the extrapolated position. Objective analyses of 
this type differ according to the scheme used in determining 
the weighting factor. One widexy used scheme, optimal inter- 
polation (01) , obtains a weight function through autocorrela- 
tion techniques (Alaka and Elvander , 1972). 01 is used 

widely for meteorological analysis (Bergman, 1979; McPherson, 
et al. , 1979) and is being adapted for ocean thermal analysis 
applications at FNOC (Innis and Williams, 1983). 
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A Navy ship at sea may not have a large number of 
observations available to generate a regular lattice of 
representative values for its operational area. Thus, util- 
izing a single observation to estimate oceanic conditions at 
distant locations is a desirable approach. The simplistic 
method of enhancing climatology described earlier is really 
a special case of 01 when only one observation is used. In 
this thesis, a method of extrapolating an observation in 
real-time to a remote position will be studied. This method 
will be referred to as "simple enhancement. " 

B. OCEANIC TIME AND SPACE SCALES 

The accuracy of ocean thermal structure estimates 
obtained using simple enhancement depends a great deal on 
homogeneity of the oceanic region under consideration. 
Knowledgeable use of the observations thus requires under- 
standing of temporal and spatial scales of oceanic 
variability . 

Significant ocean thermal structure anomalies can range 
in size from 100 km to the size of the ocean basin (TOPEX 
Science Working Group, 1981) . For instance, in the North 
Pacific, large oceanic and atmospheric anomalies have hori- 
zontal dimensions the order of 1/3 to 1/2 the size of the 
basin (Namias, 1972). More recent studies by White and 
Bernstein (1979) used autocorrelation analysis to determine 
the space and time scales of variability in the North Pacific. 
Zonal length scales were found to be 1500 km with the 
meridional scales approximately half that value. The zonal 
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length scales decrease to about 300 km at 300 m depth in the 
western portion of the mid-latitude North Pacific. This 
indicates greater variability in the western region. 
Meridional wave number spectra of temperature from hydro- 
graphic sections in this region (White and Meyers, 1982) show 
similar results of 400 to 750 km length scales at depths of 
100 to 300 m. 

These statistical results give only a crude indication 
of the horizontal space scales for which simple enhancement 
may be useful in ASW applications, i.e., several hundred 
kilometers or more. The existence of mesoscale eddies in the 
North Pacific (Bernstein, 1974; and Bernstein et al , 1982) 
will complicate the method. Large errors would be expected 
when extrapolating ' an observation across the boundary of such 
an anomaly to obtain an enhanced temperature profile. Further 
sources of error are realized as a result of relatively short 
time scale changes in the vertical structure. 

Urick (1975) noted the velocity of sound near the surface 
is sensitive to local changes in the temperature profile due 
to heating, cooling and wind mixing. Watt and Morrice (1980) 
discussed the tactical significance of sea state and near- 
surface temperature and salinity structure and noted its 
dependency on atmospheric conditions. They also demonstrated 
the significant difference in sound ray paths between morning 
and afternoon under conditions of light winds and strong 
surface warming. Intense surface thermal stratification due bo 
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diurnal solar heating is also documented by Shonting (1964) , 
who found variations in sound velocity of 2 m per second in 
the upper 5 m of the water column. 

Investigations on upper ocean response to wind mixing 
(Elsberry and Raney, 1978) show variations in mixed layer 
depth (MLD) and sea surface temperature (SST) with time 
scales of several days. Similar studies by Elsberry and Camp 
(1978) show large changes in SST are normally accompanied by 
changes in MLD and occur in association with periods of 
strong atmospheric forcing. Events of these types can 
produce significant departures from climatology (Camp and 
Elsberry, 1978) . 

The feasibility of simple enhancement can now be more 
clearly examined. Based on the preceding discussion, the 
method may seem incapable of demonstrating any skill in 
accurately depicting a temperature profile. The real skill 
may become more apparent in depicting a "representative" 
temperature profile for a given region. This may be more 
appropriate for ASW applications . It is of significance that 
active sonar systems with typical ranges of 30-40 km make an 
implicit assumption that the ocean is homogeneous over these 
ranges. Thus, range- independent acoustic models usually 
require a sound velocity profile at a single point. Due to 
internal waves and other small scale fluctuations, an 
observation of a temperature profile at a single point may not. 
be representative. The question becomes, "Is an accurate 
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depiction of MLD critical or would an average MLD for a 
region be more suitable?". Here the point can be made that 
the usefulness of simple enhancement will depend on the type 
of application for which it will be used. If only a single 
observation is available, then an average MLD cannot be 
obtained. As additional observations are taken, however, 
increasingly more sophisticated methods can be utilized to 
obtain estimates of the ocean thermal structure. The experi- 
ments in this thesis will focus on the use of a single obser- 
vation. The last experiment will address the feasibility of 
obtaining a "representative" temperature profile by using a 
filtering technique described in Appendix A. 
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II. FORMULATION OF SIMPLE ENHANCEMENT 



A. USE OF A SINGLE OBSERVATION 

A single observation, Tm(l), at position 1 can be used 
to enhance the climatology, C(x) at some other location, x, 
by the following formula: 

Ce (x) - C (x) + [ Tm ( 1 ) - C (1) ] (1) 

where Ce is the enhanced climatology and Tm consists of 
observational errors, c, plus the true value, T(l). That is, 
Tm(l) = T ( 1) + e. An alternate relation is: 

Ce (x) = Tm ( 1) + SC , 6C = C(x) - C(l). (2) 

Here it is assumed that a suitable climatology is available 
and a value at any location, x, can be interpolated from the 
climatology. Also, a real-time ocean thermal structure 
analysis or forecast (Clancy and Poliak, 1983) could just as 
easily be substituted for climatology as the "trial value." 

Equation (1) provides a convenient means of enhancing 
climatology using an observed anomaly. The equivalent form 
shown in equation (2) indicates the dependence of simple 
enhancement on the trial value used. For example, compare 
the use of a synoptic forecast with that of climatology as 
trial values. Two cases can be examined. Case (1) : The 

forecast and climatology fields may be dissimilar in magnitude 
but show the same general trends of large horizontal gradients 
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That is, one field may show a nearly consistent bias over the 
other. In this case, the SC term would be similar for both 
the climatology and forecast. Thus the values for both 
would show little difference. Case (2) : The forecast field 

is represented on a grid finer than the climatology, and 
provides greater detail than the coarser grid. The SC term 
would then be different for both fields, and the enhanced 
climatology and enhanced forecast would be different. 

B. ESTIMATING THE ERRORS USING OPTIMUM INTERPOLATION 

The extrapolation error, I , which results from simple 

X 

enhancement is: 

I = T(x) - Ce (x) , (3) 

where T(x) is the actual temperature. This relation will 
prove useful in experiments testing the feasibility of simple 
enhancement? however, a means of theoretically estimating 
the expected errors prior to performing experiments would be 
desirable. Thus the theory of optimum interpolation which 
minimizes the mean square interpolation error is examined. 

The basic equations (Alaka and Elvander, 1972) are: 

n 

T (x)_ - C (x) = S { [Tm(j) -C( j) ] P (j)} + Io (4) 

j = l 

where Io is the interpolation error and n the number of 
observations. The weights, P(j) , give the relative importance 
of each measurement, and are determined by: 
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( 5 ) 



2 

Z P(j) y (i,j) + X (i)P(i) = y (x,i) i = 1,2. ..n 

j = l 

and the minimum mean square error, E, is: 

2 n 

E = a [1 - Z y(x,i) P(i) ] (6) 

i=l 

2 

where o is the signal variance, 1/X is the signal to noise 
ratio, and y(x,i) is the autocorrelation between T(x) and 
T(i). Finally, y(i,j) is the autocorrelation between T(i) 
and T(j) , which involves products of all pairs of observations 
For n=l , equations (4) -(6) reduce to: 

T (x) = C (x) + [Tm ( 1) - C(l)] P + Io (7) 

P = y (x, 1) / (1+X 2 ) , where y(l,l) = 1 (8) 

and 

E = a 2 [l - y (x , 1) P] (9) 

with subscripts for P and X dropped. Setting P = 1 in (7) and 

comparing to (1) shows that simple enhancement is a special 

case of optimum interpolation. From (8) , P = 1 implies 
2 

y(x,l)=l and X =0. The assumption is that the correlation 
between point 1 and point x is perfect and the error of the 
observation is zero. For simple enhancement, 01 provides no 
means of estimating the errors since (9) reduces to the 
trivial solution of zero. If a value other than one is calcu- 
lated for P, the method can no longer be called simple enhance 
ment. This variation of simple enhancement with P =1 and n = 1 
will be called "01 enhancement." 
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III. 



EXPERIMENTS 



A. GENERAL APPROACH 

The concept of simple enhancement was tested by using two 
different sources of synoptic fields for trial values. 
Expendable bathythermograph ( XBT) data were used to provide 
"locally observed" temperature profiles to be extrapolated to 
the enhanced position. The data set also provided a control 
for comparisons against enhanced values. Equation (1) was 
applied and the results evaluated with equation (3). To test 
dependence on the trial value, both a climatology and a synop- 
tic ocean thermal forecast were used in the experiments. 

1. The Data 

The data used in this study were provided by the Naval 
Oceanographic Office. They consisted of a subset of XBT drops 
made by the USNS Silas Bent while surveying from Kodiak, 

Alaska to Hawaii. The sampling started on 26 September 1982 
at 54°N, 149° 30' W, continued on a track south to 
41° 45'N, and ended 1 October 1982. Observations were made 
approximately every 15 km using 750 m Sippican XBTs with 
temperature and depth accuracies of 0.2°C and 1 percent, 
respectively. Data were recorded on magnetic tape and analog 
recorders . 

Initial reduction of the data was done at the Naval 
Oceanographic Office giving temperatures at 1 m depth inter- 
vals. For this study, additional screening of the data was 



16 



made by plotting temperature profiles of the 85 XBT drops 
made. Only one XBT observation was removed from the set based 
on what appeared to be "wire stretch. " This is a condition 
in which the unravelling wire of the instrument is stretched, 
causing increased resistance in the wire and spuriously high 
temperatures . 

A grid was constructed that divided the north-south 
track into 100 equal spaces and the data were linearly inter- 
polated to it. The interpolated reports are numbered 1-100 
from north to south for reference purposes. 

The resulting vertical temperature cross-section from 
the surface to 300 m is shown in Figure 1. Large horizontal 
and vertical variability is apparent in the upper portion of 
the seasonal thermocline near 50 m. This variability is 
presumably due to shallow mesoscale eddies , fronts and inter- 
nal waves. The domain includes most of the subarctic transi- 
tion zone which has been studied extensively by Roden (1970, 
1971, 1977) . The colder water below 100 m is typical of the 
intermediate water mass of the Pacific (Reid, 1965). A 
common characteristic of the data is the appearance of temper- 
ature inversions between 100 m and 150 m (Figure 1) . These 
inversions agree with observations by Roden (1977) in nearly 
the same region. 

The location of the subarctic front is fairly persist- 
ent on a weekly time scale (Roden, 1977) . Typical horizontal 
temperature gradients in the frontal region are 1-2°C/100 km. 
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LnTITUDE 

Figure 1. Meridional temperature transect along longitude 149° 30' W. XBT 
observations made from USNS Silas Bent from 26 September to 1 October. 
Contour interval is 0.5°C. Numbers and tic marks at top show position 
of interpolated observations. Contour labels show temperature times ten. 



Two such fronts with gradients of 1.9°C/100 km occur between 
about 45° 20' N to 45° 45' N and between 43° 10' to 43° 45' N 

(Figure 1) . 

Above and below the thermocline there is almost no 
vertical coherence in the thermal field. Except for the 
variations in the 40 m and 60 m region, the horizontal length 
scale is not easily determined by inspection. However, the 
scale below 60 m can be estimated to be at least the length 
of the survey track (about 1200 km). This is consistent with 
length scales in the eastern North Pacific described earlier. 
In all respects, the data set is typical for this season and 
region of the ocean and should be sufficient for a feasi- 
bility test of simple enhancement. 

2. The Thermal Fields 

The thermal fields used to provide trial values in 
this study were a monthly climatology and a daily forecast 
of th_e northern hemisphere upper ocean thermal structure. 

Both fields are represented on the northern hemisphere 63 x 
63 polar stereographic projection used at FNOC. 

Climatological temperature profiles were extracted 
from the monthly climatology used at FNOC. The climatology 
was constructed from the surface to 400 m using data from 
the National Oceanographic Data Center (NODC) , and data 
received through direct exchange from foreign Navies and 
foreign and domestic institutions (Bauer, 1982). Its struc- 
ture is represented vertically on a fixed level grid. The 
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top of the permanent seasonal thermocline and its shape are 
described separately by a set of variable levels (Table 1) . 

Forecast temperature profiles were extracted from 
fields produced by the Thermodynamic Ocean Prediction System 
(TOPS) run daily at FNOC. TOPS cycles with a real-time ocean 
thermal structure objective analysis in an analysis-forecast- 
analysis fashion (Clancy and Poliak, 1983) and was in a test 
and evaluation phase during the period of this study. The 
vertical grid is shown in Table 1. 

The TOPS forecast presently uses the Level-2 turbu- 
lence closure of Mellor and Yamada (1974) to parameterize 
the vertical diffusion of momentum, salt and heat within the 
mixed layer. The parameter of this model that controls the 
MLD evolution is the Richardson number. The base of the 
mixed layer is the depth at which the local Richardson number 
exceeds a critical value. The Level-2 closure differs from 
those used in a bulk model, such as Garwood (1977) , where 
MLD is explicitly predicted. 

Vertical temperature cross-sections from the surface 
to 300 m (Figures 2-3) are shown for the climatology and 
forecast fields. 

Monthly climatologies for September and October were 
extracted and linear interpolation was used to obtain values 
for 29 September, the center date of the data set. Clima- 
tological profiles were then interpolated to the 100 equally 
spaced data points. This interpolation was nonlinear using 
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Table 1. 



a. Vertical grid used to represent the TOPS forecast 
temperature profiles. 

b. Vertical grid used to represent the FNOC climatology. 
In addition to fixed levels, a set of floating levels 
is used to represent the location and shape of the 
top of the seasonal thermocline in the climatology. 



a. TOPS 


b. Climatology 


Depth (M) 


Depth (M) 


Floating 


* 0 . 


0 . 


Primary layer depth (PLD) 


2.5 


25. 


Temperature at PLD 


7.5 


50. 


Temperature 25M above PLD 


12.5 


75. 


" " 12.5M below PLD 


17.5 


100. 


25 .M below PLD 


25 . 


125 . 


" " 50. M below PLD 


32.5 


150. 




40. 


200. 




50. 


250. 




62.5 


300. 




75 . 






125 . 






150. 






200 . 






300. 







*Not a computational level in TOPS numerical integration. 
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a standard subroutine (INTRPS) residing in an FNOC subroutine 
library. Finally, linear interpolation was used vertically 
to obtain values at 1 m intervals to match the data. The 
climatological profiles were then subtracted from the observed 
profiles to obtain anomaly profiles. A similar procedure was 
used to compute anomaly profiles for the TOPS forecast. 

Since the effect of diurnal and other short time scale changes 
will not be studied, a single TOPS 24-h forecast will be used 
for the entire 7 day period. During the cruise, XBT observa- 
tions were sent to FNOC and assimilated into the daily 
analysis/forecast system. Comparisons with data are nonthe- 
less independent since the TOPS forecast is valid on 26 
September, the start of the data set. 

Figures 4 and 5 show contours of the climatological 
and forecast anomalies. Note that the observed temperature 
field differs from climatology by several degrees C. Such 
anomalies are not uncommon. In the forecast, the differences 
can be explained by imperfect initial conditions, inaccurate 
atmospheric forcing, or inadequate parameterization in the 
mixed layer forecast. Regardless of the reasons for these 
deviations, their existence suggests that a technique to 
reduce the size of the anomalies, such as simple enhancement, 
would prove useful. 
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Figure 2. Same as Figure 1 but for Climatology. 
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Figure 3. Same as Figure 1 but for TOPS 24 hour forecast. 
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climatology in the three 
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Figure 5. Same as Figure 4 but for TOPS forecast anomaly. 



B. ENHANCEMENT USING A SINGLE OBSERVATION 



The 100 data points were subdivided into three over- 
lapping sections for detailed study ; points 1-40; 31-70; and 
61-100. To test simple enhancement, the anomaly profile 
(observed minus climatology) at the first point of each 
section was added to all the climatological profiles in that 
section. The first profile for each section exactly repro- 
duces the observation and adjacent altered profiles are called 
"enhanced climatology." A similar procedure was used to pro- 
duce "enhanced TOPS," except that the forecast anomaly 
(observed minus forecast) was added to all forecast profiles. 

Error fields were computed (observed minus enhanced) for 
enhanced climatology and enhanced TOPS and contoured from 
the surface to 100 m (Figures 6-11) . Large differences 
occurred between 30 m and 60 m. This indicates that the 
vertical temperature gradient at the base of the mixed layer 
cannot be accurately extrapolated forward in space using this 
method. The fact that enhanced temperature profiles for 
climatology and TOPS are nearly the same is more readily 
apparent in examination of individual profiles, which are 
shown for every fifth point in Figures 12-14. The similari- 
ties of enhanced climatology and enhanced TOPS is consistent 
with equation (2) and the discussion of Case 1 in Section II. 

Notice also that erroneous temperature inversions were 
introduced in Section 1 at points 6, 11, 16, 21, and 26 
(Figure 12). The reason for this can be explained by an 
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Figure 7. Same as Figure 6 but for points 31-70 and section 
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Figure 8. Same as Figure 6 but for points 61-100 and section 
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Figure 9. Same as Figure 6 but for enhanced TOPS forecast. 
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Figure IQ, Same as Figure 7 but for enhanced TOPS forecast. 
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Figure 11. Same as Figure 8 but for enhanced TOPS forecast. 




a sz as sa aai szt asi sat wz ^a s z 0s s l 001 szi 0s 1 sat 000, 




(U)Hid3Q (U)Hld30 





03 










CD 










O 




-P 






G 




C 






03 


• 


CD 






x: 




£ 




m 


c 


iH 


CD 




CM 


CD 


CD 


O 




0 




G 


G 


• 


Eh 


M-l 


03 


03 






0 


Qj JZ H 


03 






G 


1 


CD 


CO 


P 


CD CM 


O 


CD 


CD 




H 


c 


rH 


£ 


CD 




03 


• rH 


0 rH 


CO 


JZ 


M-l 


1 — l 


a <u 


G 


0 


— " 


£ 


p 


CD 


p 




♦H 


p 


G 


a w 


CO 


CP 


P 




c 




•H 




CD 


0 


P 


P4 


03 


P 


•H 


0 




G 


P 


-P 


M-l 


G 


03 


-P 


03 




*H 




03 


> 


-P 




>1 


P 


P 


G 


CO 


Cn 


CD 


Q) 


•H 


-P 


0 


a w 


0 


G 


1 — l 


e X3 


a 


*H 


0 


CD 


0 




0 


-P 


Eh 




CD 


a 


03 




CP O 




£ 




c 


G rH 


*H 


• 


•H 


CD rH 


»— 1 




>1 P 


03 


O 


rH 


M-l 


CD 






CD 


•H 


M-l 


M-l 


03 


G 


P 


CD 


0 


(D 


03 


CD 


P 




O 


a 


> 




G 


G 






CD 


0 


03 


u 


JC 


JC 


•H 


JZ 


CD 


-P 


-P 


-P 


c 


a 


•H 




03 


(D 


a £ 


CO 


O 


G 


p 




•H 


0 


P 


' 


03 




H 






CD 1 — ! 




M-l 


CO 


P 




P 


0 


g 


03 


-P 


0 




0 


a 


C 4-1 


CO 


■H 


£ 


♦H 




CD 


-P 


0 


0 H 


1 — 1 


03 


0 


CM 




■H 


> 






CD 


m 


P 


in 


03 


P 


0 


CD 


CM 


G 


P 


P 


CO 


0 


03 




a jq 


Eh 




•H 




0 




H 


pL| 


(U 




03 






P 


0^ CD 


G 


CD 


p 


G 


O 


0 


CD 


-p 


•H 


G 


•H 


in 


03 


>1 03 


-P 


— ' 


P 


M-l 


JZ 


O 




CD 


•H 


G 


CD 




a . 


P 


CD 


CO 


• 


£ 


CD 






G 


CD 


> 


03 


p 


0 


Eh 




G 


0 


•H 




•C 


03 


M-l 


-P 




-P 






O 


• 


•H 


>1 CD 


CD 


03 


£ 


CP 


P 


CO 


CM 




0 


03 




r— 1 


TS 


H 




CO 




CD 


0 


CO 


•H 


CD 


P 


-p 


-P JZ 


P 


03 


03 


«H 


-P 


P 


a 


£ 


P 




O' 


£ 


■H 


CO 


P 


•H 


0 rH 


CD 


0 




O 


U 


0 C 4-i 



34 




35 






CD CO 
X -H 
-p x: 
-p 
co 

•H P 
O 

( — I Mh 
ro 

-P 
■P C 
C CD 
•H g 
O CD 
CH O 
C 
03 

• x 

og c 
<D 
C 

o a> 

•H rH 

-p a 

o g 

(D -H 

co m 

nj P 

a o 

03 M-) 
•H -P 

c 

•H 

- O 

a 

ro 

<D 

- a 

rH c 



ro 


CD 




P 


CO 


CD 


-P 


MH 


C 


CD 


•H 


P 


0 




a 


<D 




x: 


p 


4J 


0 




Mn 


03 




C 


•p 


03 


3 




X 


C 




0 


05 


•H 


CM 


-P 


rH 


a 




CD 


CO 


CO 


03 


CD 


CD 


x: 


g 


-p 


fd 




w 


a 




•H 


• 


-P 


03 


c 


ro 


•H 


rH 


0 




Ou 


0) 




P 


-P ■ 


3 


CO 


tr> 


p 


■H 


-H 



a 

(D 

p p w 



36 







a sz as sz aai szi asi szi aaz a sz as sz eai szr asi szi aaz 
(U)HJLd3a C^lHidjQ 



<N 

c 

o 

•H 

■U 

O 

CD 

W 

H3 

C 

03 

LD 



m 



ID 

U) 

-p 

•H 

O 

a 

p 

c 

M-l 

-P 

.Q 

03 

CN 

i— 1 

0) 

03 

CD 

£ 

03 

01 



-Q 
no 
• — l 

CD 

U 

P 

O' 

■H 

fa 



37 




t ; 1 1 1 1 1 1 1 r 





cw)Hid3a tuwicoa 



38 



TEMPERR1URE DEG l' TEMPERED URE EJEG C TEMPERRTIIItt DEI, i 

Figure 14a. Same as 12a but for points 61, 66, 71 and section 3. Point 61 is 
the first point in the section and the reference point for simple enhancement 
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example showing the steps of simple enhancement. First note 
the depth of the inversion at point 11 (Figure 12) . A line 
indicating this depth is shown at points 1 and 11 through A, 

B, D, and E. The difference between the observation and clima- 
tology at this depth is about 2°C. Adding this difference to 
D gives the enhanced temperature E at point 11. Compare this 
to the enhanced temperature G computed the same way only 
using C, B' and F at a different depth. The apparent inversion 
due to the higher temperature at E results from the fact that 
small differences in MLD may be associated with large differ- 
ences in temperature in the MLD region. When these large 
temperature differences are added to an independent clima- 
tological temperature profile at a different location, the 
possibility of an inversion results. 

For each section, the mean and RMS errors (Tables 2-4) 
were computed for depths of 0 m, 50 m, 100 m and 150 m. 
Confidence limits were calculated using a student T distribu- 
tion with n=40. In all cases, the error is defined as the 
observed value minus the trial value. To composite the 
results , averages for all three sections over depth were com- 
puted (Table 5) . The mean errors for enhanced climatology and 
enhanced TOPS were nearly an order of magnitude smaller than 
those for unenhanced climatology. The average RMS errors for 
the enhanced values were about 1/2 that for unenhanced clima- 
tology. Also, the errors in the enhanced climatology are 
about the same magnitude as errors in the enhanced forecast. 
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Thus, two dissimilar trial values produce similar enhanced 
values with nearly the same improvements. This is a desirable 
result since the choice of the trial value used was not 
critical . 

Since simple enhancement was shown to be a special case 
of optimum interpolation, the former might be improved by 
including the known spatial statistics of ocean thermal vari- 
ability via this technique. To test this hypothesis, equations 
(7) and (8) were used to compute an 01 enhanced climatology. 

The signal to noise ratio, 1/A, and the autocorrelation, 
y(x,l) , could be computed from the 84 observations. A more 
desirable approach would be to obtain the constants from an 
independent and much larger data set. From White and Bernstein 
(19 79)., A^ = .32/. 74 = .432. The autocorrelation was deter- 
mined from 



V (x. 1) 



-c ( AX) 2 
e 



( 10 ) 



(See Clancy, 1983), where AX is the latitudinal separation 
between the observed and enhanced value. The constant c was 
determined by fitting (10) to the autocorrelation curve in 
White and Bernstein (1979). Although these statistics are 
depth dependent, values for the sea surface were used and 
applied to all depths for simplicity. Table 6 lists the 
resulting P's as a function of AX. Equation (9) was used to 
compute the 'expected minimum error for each enhanced profile. 
The average minimum error for all profiles (Table 5) is 0.62°C, 
compared to the actual mean error of 0.80°C. 
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Table 2. Mean and RMS temperature error (observation minus 
test value) for test values of climatology, enhanced clima- 
tology, enhanced TOPS and 01 enhanced climatology at depths 
of 0, 50, 100 and 150 m for section 1 (points 1-40) 



Mean Error °C With RMS Error °C With 
95% Confidence 95% Confidence 

Intervals Intervals 



DEPTH : 0 m 






CLIMATOLOGY 


•53(.43, .63) 


. 62 ( .52 , .72) 


ENHANCED CLIMATOLOGY 


-.42 (-.52,-. 32) 


. 52 ( . 42, . 62) 


ENHANCED TOPS 


-.39(-.49,-.29) 


.50 (.40, .60) 


01 ENHANCED CLIM. 


.01 (-. 08, .09) 


.26 (.18, . 35) 



DEPTH: 50 m 






CLIMATOLOGY 


- .83(-1.37,-.28) 


1.90(1.35,2.44) 


ENHANCED CLIMATOLOGY 


-1.47 (-2.02,-. 93) 


2.25(1.71,2.80) 


ENHANCED TOPS 


-1.5 8 (-2. 11, -1.04) 


2.30(1.76,2.83) 


01 ENHANCED CLIM. 


-1.1 8 (-1.73, -.63) 


2.09(1.54,2.64) 



DEPTH: 100 m 






CLIMATOLOGY 


-.71 (-.76, -.66) 


.73(.68, .78) 


ENHANCED CLIMATOLOGY 


1 

• 

o 

00 

T 

o 

00 

• 

o 

u> 


.17 (.11, .22) 


ENHANCED TOPS 


.10 (.05, .14) 


. 17 ( . 12 , . 21) 


01 ENHANCED CLIM. 


-.33(-.40,-.26) 


. 40 ( . 33, .47) 



DEPTH: 150 m 






CLIMATOLOGY 


-.54 (-. 59, -.49) 


.56 ( .51, .61) 


ENHANCED CLIMATOLOGY 


.18 ( . 13, .22) 


.23( .18, .28) 


ENHANCED TOPS 


.33 ( . 30, .37) 


.35 (.32,. 39) 


01 ENHANCED CLIM. 


-.14 (-.21, -.07) 


.26 (.19, . 33) 
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Table 3. Same as Table 2 except for points 31-70 and 
section 2 . 



Mean Error °C With RMS Error °C With 
95% Confidence 95% Confidence 

Intervals Intervals 



DEPTH : 0 m 






CLIMATOLOGY 


. 44 ( . 40, . 48) 


.46 (.42,. 50) 


ENHANCED CLIMATOLOGY 


.11 (.07,. 15) 


.17 (.12, .20) 


ENHANCED TOPS 


.19(.14, .24) 


.25 (.20, . 31) 


01 ENHANCED CLIM. 


.26 (.22, .29) 


.28(.24, .31) 



DEPTH : 50 m 






CLIMATOLOGY 


-1.11 (-1.66, -.57) 


2.04(1.49,2.59) 


ENHANCED CLIMATOLOGY 


.51 (-.04,1. 06) 


1.78(1.24,2.33) 


ENHANCED TOPS 


. 49 (-.04,1. 02) 


1.74(1.20,2.27) 


OI ENHANCED CLIM. 


- .21 (— . 80, . 37) 


1.84(1.26,2.43) 



DEPTH: 100 m 






CLIMATOLOGY 


-1.03(-1. 15,-. 90) 


1.10 (. 97,1.22) 


ENHANCED CLIMATOLOGY 


- .07(-.20, .05) 


.40 (.27, .52) 


ENHANCED TOPS 


- . 02 (-.13, . 10) 


.36 ( .24, .47) 


OI ENHANCED CLIM. 


- . 50 (— .65,— .34) 


.69 (.54, .85) 



DEPTH: 150 m 






CLimTOIOGY 


-1.10(-1 .25,- . 94) 


1.20(1.04,1.36) 


ENHANCED CLIMATOLOGY 


- .26 (-.42, -.10) 


.56 ( .40, . 72) 


ENHANCED TOPS 


- . 14 (-.28, . 01) . 


. 48 ( . 33, .63) 


OI ENHANCED CLIM. 


- .64 (-.82, -.45) 


.86(.68,1.05) 
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Table 4. Same as Table 2 except for points 61-100 and 
section 3. . 



Mean Error °C With RMS Error °C With 
95% Confidence 95% Confidence 

Intervals Intervals 



DEPTH: 0 m 






CLIMATOLOGY 


.56 (.48, .64) 


.61 ( .53, .69) 


ENHANCED CLIMATOLOGY 


.34 (.26, .41) 


.42 (.34 , .49) 


ENHANCED TOPS 


-.19 (-.29, -.09) 


.36 (.26, .46) 


01 ENHANCED CLIM. 


• .43 (.35,. 52) 


.51 ( . 42 , . 59) 



DEPTH: 50 m 






CLIMATOLOGY 


-2. 3 5 (-2. 55,-2. 15) 


2.43(2.23,2.63) 


ENHANCED CLIMATOLOGY 


.15 (-.05, .35) 


.63( .44, . 83) 


ENHANCED TOPS 


•14(-.06, .34) 


.63 ( .44 , . 83) 


OI ENHANCED CLIM. 


- .96 (-1.17, -.75) 


1.17 ( . 95,1. 39) 



DEPTH: 100 m 






CLIMATOLOGY 


-1.26 (-1. 37,-1.15) 


1.30(1.20,1.41) 


ENHANCED CLIMATOLOGY 


.38 ( .27, . 49) 


•51(.40, .62) 


ENHANCED TOPS 


.30 (.20, .41) 


. 4 4 ( . 34 , .55) 


OI ENHANCED CLIM. 


- .35 (-.46 ,-.24) 


.48(.38, .59) 



DEPTH: 150 m 






CLIMATOLOGY 


-1.41 (-1.52, -1.30) 


1.45(1.34,1.56) 


ENHANCED CLIMATOLOGY 


- .11 (-.22 , .00) 


.36 ( .25, .47) 


ENHANCED TOPS 


00 

o 

1 

OA 

(N 

• 

OA 

i — 1 
• 

1 


.38 ( .27, .48) 


OI ENHANCED CLIM. 


- .69(-.79,-.59) 


.76 ( .66, .86) 
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Table 5. Average mean and RMS error for all three sections 
and depths. 



Mean Error RMS Error 





(°C) 


(°C) 


/E 


CLIMATOLOGY 


-.596 


1.20 




ENHANCED CLIMATOLOGY 


-.057 


.67 




ENHANCED TOPS 


-.080 


.66 




01 ENHANCED CLIM. 


-.360 


00 

o 


. 62 


FILTERED ENHANCED CLIM. 


-.210 


.70 
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Table 6. Weights, P, as a function of latitudinal 
separation, AX, used to compute an 01 enhanced climatology. 



AX (Deg) 


P 


AX (Deg) 


P 


0.0 


0.6981 


2.1 


0.6040 


i — 1 
• 

o 


0.6979 


2.2 


0.5955 


0.2 


0.6972 


2.3 


0.5868 


0,3 


0.6961 


2.4 


0.5778 


0.4 


0.6945 


2.5 


0.5686 


0.5 


0.6924 


2.6 


0.5591 


0.6 


0.6899 


2.7 


0.5495 


0.7 


0.6870 


2.8 


0.5396 


0.8 


0.6836 


2.9 


0.5296 


0.9 


0.6798 


3.0 


0.5194 


1.0 


0.6756 


3.1 


0.5091 


1.1 


0.6709 


3.2 


0.4987 


1.2 


0.6659 


3.3 


0.4882 


1.3 


0.6604 


3.4 


0.4776 


1.4 


0.6546 


3.5 


0.4669 


1.5 


0.6484 


3.6 


0.4561 


1.6 


0.6418 


3.7 


0.4453 


1.7 


0.6349 


3.8 


0.4345 


1.8 


0.6276 


3.9 


0.4236 


1.9 


0.6201 


4.0 


0.4128 


2.0 


0.6122 
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The mean and RMS errors (Tables 2-5) indicate that no 
improvement is realized by including spatial statistics in 
simple enhancement. This is probably due to the use of a 
single observation for enhancement. Although the expected 
error is reasonably close to the actual error, a greater 
improvement might result if multiple observations were used 
in combination with equations (4) and (5). In addition, it 
should be appreciated that the 01 technique guarantees 
minimal errors only in an ensemble-mean sense. Thus, one 
particular realization, such as done here, might not neces- 
sarily show the advantage afforded by this technique. 

C. OBTAINING A REPRESENTATIVE ESTIMATE OF THE OCEAN 

THERMAL STRUCTURE 

The fluctuations at the base of the mixed layer noted 
earlier can be removed by a low-pass filter while retaining 
large scale trends of a few hundred kilometers or more. A 
simple procedure to "correct" each temperature profile 
using a 5-point running mean horizontally (see Appendix A) 
was applied to the data for the three sections separately. 
Anomalies from the filtered data were then used to compute 
enhanced profiles for climatology as before. Since the 
mixed layer depth was averaged over 5 adjacent temperature 
profiles, small scale fluctuations will not be present in the 
filtered data. The resulting enhanced temperature profiles 
would then be "representative" over the 5 points, or about 
48 km. The contoured error fields of enhanced climatology 
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resulting from this procedure (Figures 15-17) showed almost 
no change compared to enhanced climatology computed from 
unfiltered data. The mean and RMS errors (Tables 5, 7) also 
show no significant reduction. The results therefore do not 
indicate that filtering the data will reduce the errors of 
simple enhancement. Filtering, however, will produce 
temperature profiles which represent a larger domain than a 
profile at a single point. Although the problem of extrapo- 
lating MLD to a new location is still unresolved, the 
concept of a "representative" MLD might preclude the require- 
ment of a MLD at a single point. Testing of this idea, 
however, will not be addressed in this thesis. 

The results of this and preceding experiments suggest 
that it is feasible to extrapolate the observed temperatures 
of most of the water column to a new location. The problem 
of determining the MLD should be dealt with separately. 

Since these experiments are relevant only to a specific time 
and domain, additional studies under different conditions 
are necessary to assess properly the usefulness of simple 
enhancement. 
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Table 7. Mean and RMS temperature error of filtered 
observations minus enhanced climatology computed from 
filtered observations (see text for explanation) . 



Mean Error °C With RMS Error °C With 
95% Confidence 95% Confidence 

Intervals Intervals 



DEPTH: 0 m 






POINTS 3-38 


-.67 (-.77, -.56) 


.74 ( .63, .84) 


POINTS 33-68 


1 

• 

0 
Ln 

M 

O 

1 

• 

o 


.14(.10, .18) 


POINTS 61-96 


.32 (.24,. 41) 


.41(.32, .49) 



DEPTH: 50 m 






POINTS 3-38 


-1.67 (-2.17,-1.17) 


2.23(1.73,2.73) 


POINTS 33-68 


.19 (-. 42, . 80) 


1.82(1.21,2.43) 


POINTS 61-96 


.27 (.11, .44) 


.55 ( . 39, . 72) 



DEPTH: 100 m 






POINTS 3-38 


-.17 (-.22, -.11) 


.24 ( .18, . 30) 


POINTS 33-68 


-.41 (-.54 ,-.27) 


.56 (.43, .70) 


POINTS 61-96 


. 36 ( .25, .48) 


.51 (.39, .63) 



DEPTH: 150 m 






POINTS 3-38 


-.07 (-.12 ,-.02) 


.16 ( .11, .21) 


POINTS 33-68 


-.42 (-.59, -.25) 


CM 

00 

• 

00 

LO 

to 

• 


POINTS 61-96 


-.15 ( .26 ,-.03) 


.37 (.26, .49) 
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Figure 16. Same as Figure 15 but for points 33-68 and section 
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Figure 17. Same as Figure 15 but for points 61-96 and section 



IV. CONCLUSIONS 



A method of extrapolating an observed temperature profile 
from one location to another was investigated. The technique, 
referred to as simple enhancement, was examined as a possible 
application to real-time Navy ASW operations . The extrapola- 
tion required the use of a trial value which was obtained 
from two sources: an ocean thermal climatology and a real- 

time ocean forecast model (TOPS) . An enhanced climatological 
temperature profile was obtained by adding an observed anomaly 
(jL.e., observation minus climatology) to the climatology at 
some other desired location. The same method was used to 
obtain an enhanced forecast temperature profile. This proce- 
dure was shown to be a special case of optimum interpolation 
when only one observation is used, the noise- to-signal ratio 
is zero, and the autocorrelation function between the two 
points in question (i.e., the point of the observation and the 
point where the enhanced profile is desired) is one. 

The feasibility of simple enhancement was evaluated by 
calculating mean and RMS errors of enhanced and unenhanced 
error fields. Vertical transects of the error fields were 
contoured to show the magnitude and the distribution of errors 
associated with simple enhancement. Vertical profiles were 
also examined. 
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The results were composited by averaging the mean and 
RMS errors for each section and depth. The average mean 
error for enhanced climatological and TOPS errors was nearly 
an order of magnitude smaller than the average for unenhanced 
climatology. The average RMS errors for the enhanced fields 
were about 1/2 of that for unenhanced climatology. 

Examination of the contoured error fields showed the 
magnitude of error to be 1 to 2°C at the base of the mixed 
layer. In addition, temperature inversions were introduced 
in some of the enhanced temperature profiles. This is an 
indication that the MLD could not be properly extrapolated. 

The choice of the trial value did not appear critical. 
Enhanced vertical temperature profiles of climatology and 
TOPS appeared nearly identical , even though the unenhanced 
values were different. Simple enhancement was also tested 
using known spatial statistics to compute an observational 
weight from the Optimum Interpolation formalism. The errors 
in the resulting 01 enhanced climatology did not show a 
reduction compared to simple enhanced fields . This is 
probably explained by the use of a single observation to 
compute the enhanced temperature profile. The 01 technique 
guarantees minimum error only in an ensemble-mean sense using 
many observations. Thus, the approach of incorporating spatial 
statistics might not be appropriate when a single observation 
is used. 
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Finally, filtering the data prior to enhancement provides 
a method of producing a more representative temperature 
profile for a region. Enhanced climatology produced from 
filtered data shows no reduction in error over enhanced clima- 
tology produced from unfiltered data, but shows less error than 
unenhanced climatology. Thus, it is possible to produce an 
enhanced temperature profile that is also representative of a 
region, rather than a single point. The disadvantage of 
filtering is that more than a single observation is required. 

This method could provide the Navy with an improved ASW 
capability. Screening of the results, however, is recommended 
and suggests the development of an objective technique for 
removing erroneous temperature inversions. Variations of 
simple enhancement are possible. In the present scheme, 
temperatures at each depth are extrapolated forward in space. 

A temperature profile, however, can be defined by parameters 
other than temperature and depth pairs. A reasonable varia- 
tion would be to define a temperature profile by SST, MLD, 
and temperature gradients below the MLD. More parameters 
could be defined and added if further detail in the tempera- 
ture profile is desired. Then the observed SST, MLD and 
temperature gradients could be used to compute anomalies of 
these parameters and added to trial values of the same param- 
eters at a new location. An enhanced temperature profile 
reconstructed from these new parameters would not show temper- 
ature inversions. 
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An adequate climatology is all that is necessary for 
simple enhancement. The use of TOPS forecast profiles works 
equally well, but a better use of TOPS might be to predict 
changes in MLD and SST. This would reduce errors in short 
time scale changes by knowing if conditions will persist or 
change in time. 
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APPENDIX A 



FILTERING THE DATA 

Consider a small scale fluctuation at depth Z (Figure 18) 
displaced from the mean depth, Z , by a distance DZ m . Iso- 
therms at depth ZO above and below Z oscillate about their 
own mean positions Z0 m , with departures, DZg, that decay 
toward zero at the surface and some depth below Z. The 
effects of the fluctuations can be removed by finding Z Q , the 
depth of each isotherm and replacing it with the mean depth, 
ZC) , of that isotherm. Assume the departures from Z are 
proportional to departures from Z0 m - Then DZ m and DZg are 
related by : 

DZ Q = F(Z)DZ m (11) 

An appropriate F(Z) is: 



F(Z) 




Zn>Z 

0 m 


(12) 


F ( Z) 


= z o/ z m 


Zg < z m 


(13) 



where D is the depth to which the effects decay to zero. 

Notice that F(Z) goes to zero as Zg approaches D, or as 
Zg approaches the surface. By definition, 

DZg — Z0 m — Zg , and (14) 

DZ = Z - Z. (15) 

m m 
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Figure 18. Sinusoidal representation of small scale 
fluctuations at base of mixed layer. (See text for 
explanation of variables.) 
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Then : 



Z °m * z 0 + <z 0/ 3 m> < z m - z > 




(16) 



Zo m - Zq + (D-Zq) ( z m — Z ) / (D Z m ) Zq > z m 



(17) 



This formal approach is easily applied. First, select an iso- 
therm at the base of the mixed layer. Determine the mean 
depth, Z , over five adjacent temperature profiles using a 
five point running mean horizontally. Next find Z, the depth 
of the selected isotherm for each profile. Then for each 
level, Zq, compute the corrected depth Zq^ using equation (16) 
or (17). The resulting profile will now have new corrected 
depths at each level. The temperatures remain unchanged and 
only their vertical position has been altered. 
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